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Abstract 

A system of a metastable phase with several sorts of the heterogeneous centers is considered. An analytical 
theory for the process of decay in such a system has been constructed. The free energy of formation of the 
critical embryo is assumed to be known in the macroscopic approach as well as the energy of solvatation. 
The process is split into some periods and the analytical description of every period is given. The most 
difficult to describe is the period of the essential formation of the embryos of a new phase. At first some 
asymptotes are investigated and then a general solution is suggested. Several approximate transformations 
are accomplished with the corresponding numerical estimates and some analytical justifications. 
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1 Introduction 



The case of condensation of the supersaturated vapor into the state of liquid droplets seems to be the simples 
one among the first order phase transitions. Traditionally this case is treated as a model to introduce some nei 
theoretical constructions in the description of the first order phase transitions. This leading role of the cas 
of condensation was outlined by the creation of the classical theory of nucleation by Volmer 0, Becker an 
Doering @, Zeldovitch g and Frenkel §J which gave for the first time a rather simple expression for the rat 
of nucleation (i.e. for the rate of the appearance of new droplets). 

The careful analysis of the classical theory leads to the great number of the publications with the variou 
reconsiderations of the classical expressions. Among them one has to notice the account of the internal degrees c 
freedom in the embryo made by Lothe and Pound ||, another modifications made by Reiss, Cohen and Katz |^ 
Reiss (Q , Fisher M . The application of the density functional theory to the first order phase transition made b 
R.Evans and D.Oxtoby Q and D.Oxtoby and D.Zeng [|l3| allowed to put the microscopic (mesoscopic) model 
for the condensating substance as the base for the expression for the free energy of the critical embryo and fo 
the rate of nucleation. The latest reconsideration of the classical theory can be associated with the contribution 
of Reiss, Tabazadech and Talbot p5| , Reiss, Ellerby and Weakliem ||] - Q, Oxtoby and Talanquer Q. I 
the cited publications the role of the environment of the droplet is carefully analysed. The choice proposed i 
these publications is equivalent to some specific choice of the statistical ensemble. 

Nevertheless one has to notice that there is no perfect coincidence between the concrete theoretical predic 
tions0 and the experimental results. The macroscopic expression reproduces only the qualitative behavior c 
the experiment results. Meanwhile, the relative deviation between the rate of growth of the nucleation rate i: 
the theoretical predictions and in the experimental results is rather small. 

These efforts allowed to start the investigation of kinetics in the field of the first order phase transitior 
The qualitative description of the phase transition was started by Wakeshima who considered some tim 
lags for condensation. As far as the characteristic time of formation of the droplets spectrum was necessary t 
ensure the correct experimental definition of the stationary rate of nucleation namely this value was investigate 
in In (3^] the picture of the more realistic phase transition under the smooth behavior of the exterm 

conditions (i.e. under the external conditions of the dynamic type) was given. The form of the spectrum i: 
these publications was introduced by some artificial models. But the real problem is to determine the form c 
the spectrum. 

The homogeneous decay of a metastable phase was investigated in [ fl7| by Kuni, Grinin, Kabanov. Th 
analogous consideration for the heterogeneous decay was fulfilled by Kuni in Q where only two limit situation 
were considered. The consideration of the simultaneous heterogeneous and homogeneous decay was made i: 
[ fl9| , |Q by the formal generalization of the procedure for the homogeneous decay jL7| . 

The homogeneous condensation under the smooth behavior of the external conditions (i.e. under the dynami 



conditions) was investigated in , |27| . The case of the heterogeneous condensation was analyzed in 28 , p£ 
where the based qualitative descriptions of the phase transition were given. 

The external conditions of the dynamic type and the decay type are the most natural external condition 
and can be treated as some basic ones. 

Practically in the condensing system there are several types of the heterogeneous centers of a different nature 
Moreover in the process of condensation on the ions the free energy of the critical embryo depends on the sig 
of the electric charge. As it is shown in [|TJ the free energy F of the near-critical embryo in the state of th 
internal equilibrium can be presented in the following manner: 

F = -bv + c«/ 2/3 + Cl is 1/3 + (c 2 + c 3 )v- 1/3 + c In v - G (1 

Here and in the further considerations all energy-like values are expressed in the units of the thermal energy 
(kb is the Bolzman constant, T is absolute temperature); a, 6, Co, ci, C2, C3 are some constants; G is the energy c 
the nuclei solvatation. It necessary to notice that in contrast to a, 6, Co, c\, ci which don't depend on the sign c 
the charge q the value of C3 is proportional to q. The value of G also depends on the sign of q. This dependenc 
is similar to dependence (Q) of ^(without G). The only specific feature we must fulfil is to substitute the numbe 
of the molecules v e of the solvatated ion instead of v. As far as v e ^ v c for the near-critical embryo the value c 

1 Certainly, one can not check the general recipes as to calculate the objects like the statistical sums. 
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F depends on the sign of q. Hence, in the presence of the radiation we immediately have two sorts of the center 
(positive and negative) with the different heights of the activation barrier, i.e. with the different activities c 
the heterogeneous centers. 

When one has the spectrum of the sizes of the solid nucleus of condensation with the weak interaction on 
must immediately come to the spectrum of activities of the heterogeneous centers. Really, in the most simpl 
model appears when the nuclei is essentially noninteractive and one has simply to add the surface term to th 
number of the molecules imaginary contained in the volume occupied by the nuclei. The free energy of th 
critical embryo is given by 

F=-bu + a(u + ^) 2 / 3 -G (2 
6vi 

where r is the radius of the nuclei, vi is the volume occupied by a molecule of the condensated substance in th 
liquid phase. So, the spectrum of sizes of the heterogeneous centers immediately initiates the spectrum of th 
activities of the heterogeneous centers^} 

Any spectrum of the activities of the heterogeneous centers can be split into several continuous parts (ma 
be also like 6- functions) which will be considered as the "types" of the heterogeneous centers. So, one can se 
four rather natural kinetic problems: 

• The process of the decay of the metastable state on the several types of the heterogeneous centers. 

• The process of the decay of the metastabtable state on the continuous spectrum of the heterogeneou 
centers. 

• The process of condensation on the several types of the heterogeneous centers under the dynamic condi 
tions. 

• The process of condensation on the continuous spectrum of the heterogeneous centers under the dynami 
conditions. 

The set of these publications is based on the ideas proposed in |32| - jj7| 

Nevertheless the theory of the heterogeneous decay was constructed only for one type of the heterogeneou 
centers. The task to construct the kinetic theory for these situations is rather essential, ft will be completel 
fulfilled here. 

The equations of condensation are rather similar in these four situations. Unfortunately, the methods c 
their solution are absolutely different. This is caused by the specific nonlinear character of the condensatio 
equations. Moreover, the analysis of any separate situation can not be spread on another one. That's why th 
four separate publications are devoted to the kinetics in these situations. 

One has to study only the periods of the intensive formation of the droplets. But in contrast to condensatio 
on the one sort of the heterogeneous centers p9| these periods aren't very shortQ 

We shall use the following physical assumptions: 

• the thermodynamic description of the critical embryo, 

• the random homogeneous space distribution of the heterogeneous centers, 

• the free-molecular regime of the droplets growth, 

• the homogeneous external conditions for the temperature and for the pressure, 

• rather a high activation barrier^, 

• the absence of the thermal effects. 



2 Later the value of the activity will be accurately defined 

3 Under the dynamic conditions 

4 The theory without the heterogeneous activation barrier is much more simple. 
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As far as according to |32| the condensation equations in the general conditions of the condensation proces 
are analogous to the case of the free molecular consumption of the vapor we shall study namely this case. 
We assume the total number of the heterogeneous centers to be constant in time. 

As far as the most interesting characteristics of this process are the numbers of the heterogeneously forme' 
droplets of the different types we shall estimate the accuracy of the theory by noticing the error of the obtaine< 
solutions for these values. The whole process of the metastable phase decay can be split into two periods: th 
period of the essential formation of the droplets and the period of the essential consumption of a metastabl 
phase. At first we shall investigate the period of formation of the droplets. The unit volume is considerec 
The consideration of the further evolution is rather simple and can be reduced to the first order diffcrcntiE 
equation in the manner of |25| for the situation of the decay and in the manner of |2^] , for the situation c 
the dynamic conditions. We needn't to consider this period here. 

In order to present the clear descriptions some common facts will be recalled at the beginning of ever 
publication. It allows to consider the theories as the separate ones. 

Below the situation of the decay of the metastable phase on the several sorts of the heterogeneous center 
will be described. 

We assume the total number of the heterogeneous centers to be constant in time. 



2 Kinetic equation 

Suppose that there are several sorts of the heterogeneous centers. We shall mark the total number of th 
heterogeneous centers by r\tot % where i corresponds to some sort of the heterogeneous centers. The real value 
of the free heterogeneous centers which may be solvatated but aren't occupied by the super-critical embryc 
are marked by the value of r\i. The index i or j below the value marks the sort of the heterogeneous centers 
The sort of the droplets means the sort of the heterogeneous centers. The absence of this index points tha 
the formula is valid for an arbitrary sort of the heterogeneous centers. The density of the molecules in th 
equilibrium vapor is marked by rioc, the density of the molecules in the real vapor is marked by the value r, 
The power of the metastability will be characterized by the value of the supersaturation 

_ n-rioa 



We shall define the super-critical embryos as the " droplets" . Every droplet is described by the molecule numbe 
v , or by the linear size p — v 1 / 3 . Due to the free-molecular regime of the vapor consumption we have 

d P > 

where a is the condensation coefficient and r is some characteristic collision time obtained from the gas kineti 
theory. 

The characteristic time £* will be the time of the beginning of the process of condensation. The values a 
the moment i* will be marked by the lower index t*. 

The frontal type of the size spectrum allows to introduce the frontal size z according to 

z= j (ar^dt' (3 
Jo 

Until the coalescence |2^] , [^3| which isn't considered here equation (||) ensures the growth of z in time and ca: 
be inverted 

t(z)=[ Z ra-^ (4 
Jo CW 

Hence, all values dependent on time become the values dependent on z and the relative size x — z — p can b 
introduced. During the whole evolution the droplet has one and the same value of the variable x. Considerin 
the value t(x) as the moment when the droplet with the given x has been formed (as a droplet) we can conside 
the functions of time as the functions of x . Hence, we can see that the kinetic equation is reduced to the fac 
that every droplet keeps the constant value of x. To reconstruct the picture of the evolution one must establis 
the dependencies t(z) and ((x). 
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3 Condensation equations system 



We shall mark by the argument oo the total values of the characteristics formed during the whole condensatio 
process. 

Immediately after the creation the value of the supersaturation falls down to the value 

^ = C(Q) _ ZiVtotiVei (; . 
77-00 

where v e j is the number of the molecules of the condensated substance in the equilibrium heterogeneous embryc 
During the period of the essential formation of the droplets one can assume that the value v e j is the constan 
one and take it at ( — Ordinary v e j can be taken at £ = C(0)- 

The following statements can be analytically proved for the consideration of this process: 



(1) The main role in the vapor consumption during the evolution is played by the super-critical embryos 
i.e. by the droplets. 

(2) The quasistationary approximation for the nucleation rate is valid during the period of the essentk 
formation of the droplets. 



The justification of the second statement uses the estimate for the times t\ of the establishing of the stationar 
state in the near-critical region which can be found in || , |24|] (for the heterogeneous barrier the consideration i 
the same one). Here it is necessary to remark that, certainly, there may exists some rather huge times t\ . The 
correspond to some rather big values of the number of the molecules inside the critical embryos. The half-widt. 
of the near-critical region can be estimated by the homogeneous value at the supersaturation correspondin 

2/3 

to the same value of the number of the molecules inside the critical embryo. It is proportional to v ' ■ . Her 
and in the further considerations the lower index "c" marks the values for the critical embryos. As far as th 

2/3 • •• •• • 2 /3 

absorption ability is proportional to v c i and the size of the near-critical region is proportional to v c th 
value of t\ is proportional to v 2 Jf ■ The big value of the activation barrier AFi = Fi{v c ) in the case when v c 
is greater than the characteristic length Aa: of the size spectrum means that these kinds of the heterogeneou 
centers are excluded from the kinetic process. 

One can analytically prove that during the period of the intensive formation of the droplets of one sort o 
for all sorts 

• the variation of the supersaturation allows the estimate 

IC-C*I<^ 

rfA.F(C) . 

li = -<P* — — c= $„ 6 

For the majority of the types of the heterogeneous centers the following approximations of the nucleatio 
rates Ji are valid during the period of the essential formation of the droplets 

Ji = Mvtot u gjexpfo fc"**) )-^ (7 

Vtot i 

The validity of these approximations is justified for the heterogeneous embryos with the monotonous interac 
tion between the center and the molecules of the condensated substance weaker or equal than the interactio 
reciprocal to the space distance. 

Let /* i be the amplitude value of the distribution of sizes of the heterogeneously formed droplets measure- 
in units of rioo. The supersaturation and the number of the heterogeneous centers r\tot % are the ahead 
known values. Then the stationary distribution / can be easily calculated by the following known formulas: 

JiT _._ W+expj-A^r 



where 



a<>oo n 00 n 1 / 2 A e ii/A c tvQa 
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where W + is the number of the molecules absorbed by the critical embryo in the unit of time, A e v is th 
characteristic width of the equilibrium distribution 

i/=(i/ c +i/ e )/2 

and A c (u is the halfwidth of the near-critical region 

2 l/2 

A r V = 



Q2 F ,1/2 

7h7? \v=v c 



We shall mark by n^gi the total number of the vapor molecules in the droplets formed on the centers c 
the sort " i" . To simplify the formulas we shall use 



Vtot i 



Using the conservation laws for the heterogeneous centers and for the molecules of the substance we obtai 
for gi, 9i the following equations 

9l = /* i !\z - xf exp(-T i ^^)e i dx ee G.fVj^,) (9 
Jo 4>* . 

di = eM-f* f exp(-r^^)dx) ee SiC^Qj) (10 



Vtot i Jo 

J 

where /» i = Ji{r\tot i, $*)r/$*an 00 . 

These equations form the closed system of the condensation equations. This system will be the subject c 
our investigation. For simplicity we shall investigate it for i = 1,2. As far as we measure the accuracy of th 
theory in the terms of the error in the droplets number we define these values as the following ones: 

N i = r ho ti(l-0i(z)) = Qi(0i) (11 

The spectrum of sizes can be found as the following one 

fi = UieM-T^i^)Qi (12 



4 Formal generalization of iteration method 

The formal transformation of the iteration method for the process of condensation on the one sort of th 
heterogeneous centers leads to the following equations: 



9i (i+i) - ^»(y^gj (0'^» (o) ( i3 

3 

3 

N i{ i ) =Q i {6 i{l) ) (15 

9i (o) = (16 

0i (o) = 1 (17 

9i (l) = /* i-r- (18 



G 



0j (i) = exp(-/» (19 

Vtot i 

Ni (2) (oo) - r? tot 41 - cx P (-.A ^)" 1/4 ^)] (2° 

where 

/■OO 

A= / cxp(-a; 4 )dx w 0.9 (21 
Jo 

The third iteration can not be calculated in the analytical form. Meanwhile one can prove the convergence c 
the iterations analytically. 

Let us analyze the expression for AT i ( 2 )(oo). Assume for a moment that for some i and j 

/* i ^> /* j 

Let us decrease rjtot % keeping the constant value of /* j which is proportional to rji with the fixed dependenc 
on ( by the increasing of the activity of the heterogeneous centers of the sort "i", i.e. by the decreasing of th 
height of the heterogeneous activation barrier. It is obvious that when r/ tot j is small then the total numbc 
of the heterogeneously formed droplets coincides with the total quantity of the heterogeneous centers and goc 
to zero when r\ to t % goes to zero. The value of gi at the end of the period of formation of the droplets on th 
heterogeneous centers of the sort "j" can be estimated as 

yi — 

where Axj is the width of the size spectrum (the size distribution function) of the droplets of the sort "j". Th 
value of Axj is above restricted by the value Axj which is the width of the size spectrum without any influcnc 
of the droplets of the other sorts. Certainly, the value of Axj doesn't depend on rjtot % and on /„ , for all ; 
Hence, the influence of the heterogeneous centers of the sort i on the process of condensation on the centers c 
the sort j becomes unessential (negligible) in the limit rj tot i — ► 0. At the same time the expression for Nj (2)(oc 
shows that in the limit r\ tot j — > 0, /* j = const the influence of the heterogeneously formed droplets of a soi 
i doesn't become unessential. This leads to the enormous error in Nj (oo). One cannot obtain the analyticc 
expression in the third approximation for Nj in the frames of the standard iteration method and the secon 
iteration gives the wrong qualitative results. 

The reason for the deviation in the results from the pure heterogeneous consideration is the following one 
In the case when the interruption of the embryos formation is caused by the exhaustion of the heterogeneou 
centers the error in the value of g is compensated by the squeezing force of the operator 5,. The analogou 
property isn't valid for the operator Qi due to the cross influence of the droplets formed on some different sorts 

Note that the procedure of this section can be successfully applied for the heterogeneous condensation o: 
the centers of the one sort. Then one can analytically prove that the relative error of N( 2 ) is less than 0.015. 

5 Limit cases 

The remarkable property of condensation in the situation considered here is that we can construct the theor 
by the simple investigation of all limit cases and cover practically all situations. 

At first we shall extract the characteristic sizes in order to construct rather a simple solution. 

5.1 Characteristic lengths 

The direct generalization of the iteration method fails due to the unappropriated account of the cross influence c 
the heterogeneous droplets of some different sorts. Nevertheless it allows to define the spectrum of the droplet 
in a proper way when the cross influence is eliminated. Hence, we may use it to get the characteristics of th 
" self- formation" of the droplets of the various sorts. 
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On a base of the first iterations in the general procedure one can realize that for some process of a separat 
formation there are two characteristic lengths. The first one is the length of spectrum in the situation whe: 
there are no exhaustion of the heterogeneous centers (and no droplets of the other sort). One can say tha 
condensation occurs in the pseudo-homogeneous way. For this characteristic value we have: 

A^ = (^) 1/4 (22 

J- if* i 

This length is going from the first iteration for if the process is split into the separate processes of condensatior 
The second length is the length of the spectrum when there is no vapor exhaustion but only the exhaustion c 
the heterogeneous centers. Then the width of the spectrum is: 

r Vtot i ,„„ 

OiX = (23 

This length is going from the first iteration for 6^. 

Practically the hierarchy between AiX,8jX is ensured by the hierarchy between /» 4,77* j. The values of I 
are rather (in comparison with /* j) unsensible to the value of the supersaturation. Really: 

Ti _ dAF dF c dG 

In the frameworks of the barrier character of condensation (AF ^> 1) we can give the above estimate fo 
(dF c /d() for the force of the interaction between the heterogeneous center and the molecules of a liquid whe 
this force decreases monotonously in space. The value of (dF c /d() can be above estimated by it's value i 
the limit of the homogeneous condensation (dF c tom/dC,). As far as the energy of solvatation depends on th 
supersaturation essentially weaker than {dF c /dC,) we can neglect the last term of the previous equation an 
obtain: 

dAF dF c hom 

~c dT~ (25 

But this dependence is rather weak one in comparison with the very sharp /* j dependence on the supersati; 
ration. 

Another one important fact is the frontal character of the back side of the spectrum in the pseudc 
homogeneous situation (when I\ plays some role)^]. The frontal character can be seen from 

p. 

fi = /* i e M-^-(Yl f* ( 26 

* 3 

Hence, the essential variation of the length A^a; can be caused only by the enormous variation of /* j. 
Instead of SiX we shall use the parameter 

hi = ^ (27 
AiX 

In the reasons of simplicity we shall restrict ourselves by two sorts of the heterogeneous centers. 



5.2 The case A±x ~ A 2 x 
5.2.1 Situation hi < 1, h 2 > I 

In this situation the following fact can be noticed: 

• The process of formation of the droplets on the heterogeneous centers of the first sort (the "first sor 
droplets") doesn't depend on formation of the droplets of the second sort. 

5 The pseudo homogeneous situation can be treated as the situation when the centers of condensation remain practically une> 
hausted (in relative sense). 
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It can be directly seen from the chain of inequalities 

Six < Aix ~ A 2 x < S 2 x (28 
So we can describe the process of formation of the first sort droplets by the following equalities: 

9i = Ui I (z-xfcxp(-r 1 ^^)9idx = Gi(g 1 ,6i) (29 
Jo 

6 1 = exp(-/* [ Z exp (-Ti^)dx) = SiO/i) (30 

This system can be solved by the following iterations 

9l (i+l) = Gl(5! (31 

#i (»+i) = Si(gi w ) (32 
W 1(i) =Qi(0 1(<) ) (33 
The operators Gi, 5i and Qi have some remarkable properties. When for all values of the arguments we havt 

Wi < W2 

then 

Si(wi) < Si(w 2 ) 

is valid for all values of the arguments. When for all values of the arguments we have 

Wi < w 2 

then 

Qi(wi) > Qi{w 2 ) 

is valid for all values of the arguments. When for all values of the arguments we have 

Wi < w 2 
Vi > v 2 

then 

Gi(vi,wi) < Gi(v 2 ,w 2 ) 

is valid for all values of the arguments. Let us choose the initial approximations as the following ones 

9i (0) = (34 
0i (o) = 1 (35 

We see that 

9i (o) < 9i 
9i (o) < 9i (i) 

for an arbitrary i, and 

#1 (o) > 
&i (o) > ^1 (») 

for an arbitrary i. Particularly, the following estimate is valid 

Ni (2) < Ni < Ni (3) (36 
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These estimates allow to prove the convergence of the iterations. The calculation of the iterations gives 

z A 

9i (l) = /* i— (37 
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0i (i) = exp(-/* i— ^-z) (38 

Vtot 1 

N, (2) (oo) = r? tot i(l - exp(-/, i-^(-^)- 1 /4 / -V^ )) (39 



%ot i 4$* 

As far as one can prove that 



^-\N 1{i) -N 1{j) \>0 (40 



dx 

then by the simple numerical calculation of A r i( 3 )(oo) one can easy obtain that 

Ni ( 2 ) - -ZVi I 

( 3 - < 0.015 



which completes the iteration procedure. 

On the base of the iterations one can get some approximations for the supersaturation: 

C(i+i) =*•-/. i J (z-x) 3 cxp(-ri^)0i {l) dx (41 

The strong inequality allows to obtain the second approximation for ( 

C(2) =$*-/* l / (z- xf exp(-Hx)dx (42 



where 

H = (43 

Vtot 1 

or after the calculation 

z 3 3z 2 6z 6 6 
C( 2 ) = + /. + -jjt - + - jji cxp(-Hz)} (44 

This expression can be simplified. Let us notice that the supersaturation appears in the expression for th 
size spectrum f(x) in the following form 

exp(-r^) 

After the substitution of Q2) into this expression we realize that in the case when £ deviates essentially fror 
$ all the terms except the first two terms can be neglected: 

C( 2 ) = - z^-fl (45 

Hence, for the second sort we can obtain the following system of equations 

92 = J* 2 / {z - xf exp(-r 2 )9 2 dx = G 2 (g 2 + (Vtot l/n^x 6 , 2 ) 

Jo 

and 

2 = exp(-/» 2 / exp(-r 2 )dx) = S 2 (g2 + (Vtot i/rioo)^) (47 

»7ioi 2 Jo ®* 

Having introduced 

A 2 = 92 + (Vtot l/n^z 3 (48 



(46 
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we can rewrite this system as 

A 2 = G 2 {\ 2 ,6 2 ) + ( Vtot i/ noo )z 3 = G+(X 2 ,9 2 ) (49 

9 2 = S 2 {\ 2 ) (50 

The operator has the same properties as Gi,G 2 have. All estimates (|34|)-(j36|) remain valid after th 
substitution of the index "2" instead of "1" and the operator G + instead of G. Moreover one can see that 

which leads to the same estimate of the relative error (0.015) as in the pure heterogeneous condensation of 
separate sort. 

Actually we can avoid here the calculations according to such a complex procedure. Let us notice that th 
term vtot iz 3 jn^ ensures the characteristic length 

Dl = (^2S_)V3 (52 

1 2Vtot 1 

Really 

Di > eA lX ~ eA 2 x e ~ (3 ± 1) (53 

So, condensation on the centers of the second sort occurs in the separate way and we can use formulas (p9])-(|33] 
©-© witn tne substitution of the index "2" instead of the index "1". 

5.2.2 Situation h\>l, /ij<l 

As far as Ai ~ A 2 we can change the numbers of sorts and reduce this situation to the previous one. 

5.2.3 Situation hi > 1, h 2 > 1 

In order to analyze this situation we must realize why in the separate condensation of one sort the second iteratio 
gives rather precise results. This fact is explained by the free molecular regime of the vapor consumption whic 
leads to the power 3 in the expression for g. Due to rather a great (in comparison with 1) value of this powe 
the droplets of the big sizes situated near the front side of the spectrum are the main consumers of the vapoi 
These droplets have rather a small (in comparison with Aix) values of the variable x. The exhaustion of th 
heterogeneous centers due to hi > 1, h-z > 1 doesn't affect in a strong manner on the process of their formatior 
Hence, the cross influence is rather weak and we can use the general iteration procedure. After the calculatio 
of the iterations we have 

91 (l) = /* l-j- (54 
6i (n = exp(-/* i—^-z) (55 

Vtot 1 

Ni (2) (oo) = vtot i[l - e*p(-/. i — ((^)/* i + (^-)/* 2)- 1/A A)] (56 

Vtot 1 4 < J>» 

2 4 

92 (i) =1*2^ (57 
6» 2 (i) = exp(-/* 2 -^^z) (58 

Vtot 2 

N 2 (2) M = vtot 2 [1 - exp(-/, 2 ^((^)/» i + (if)/. 2)- 1/i A)} (59 
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5.2.4 Situation hi < 1, hi < 1 

As far as 

Six <C A x x - A 2 x (60 

the droplets formed on the heterogeneous centers of the second sort don't act on the pricess of formation of th 
droplets on the heterogeneous centers of the first sort. Due to 

S 2 < A 2 x - An (61 

the same fact is valid for the droplets formed on the other sort of the heterogeneous centers. Hence, th 
system is split into some parts corresponding to the separate processes of condensation on the different sorts 
Constructions (^9|)-(^3|) can be reproduced here. But in the case hi <C 1 for all i the vapor exhaustion can b 
neglected in comparison with the exhaustion of the heterogeneous centers. 

In this situation we can obtain some precise explicit results for the values of fi(x), 0i(x): 

1 (x)=exj>(-f*i—z) (62 

Vtot 1 
71 

9 2 (x) = exp(-/» 2 — —z) (63 

Vtot 2 

77 

/i(aO = /.exp(-/. i— *) (64 

Vtot 1 

/a(aO = /* exp(-/« 2 — —z) (65 

Vtot 2 
^oo Ttoc 

This expression for ( is obtained by the same procedure as that which led to (^5|). 
5.3 Case A x x < A 2 x 

The case A 2 x ^> A\x is simmetrical to this case and can be considered by the simple change of the indexes. 

Due to Aix <C A 2 x the droplets of the second sort don't act on the process of formation of the droplets c 
the first sort. Hence, the process of formation of the droplets of the first sort can be described by the iteratio 
procedure from the section 5.2.1 (eq.(p^-(f40|)). 

5.3.1 Situation hi < l,h 2 > 1 

Due to hi -c 1 the equation for the first sort can be simplified and we have the following equations 

,9i = /* i f (z- xf exp(~Hx)dx ~ 7 J^H Z 3 (67 



The value of Bi is given by (62), the value of fi is given by (|64j). 

For condensation on the centers of the second sort we have some equations analogous to (|46|) , (fl7|) . So w 
can adopt here equations (p8|)-(]5l|). But in this situation inequality ( p3| ) isn't valid and we must calculate th 
iterations. We can choose X 2 ( ) = and get 

4 

A 2(1) =/*4 + ^-V (68 
2 (2) = exp(-/. 2 -^L_ f exp(-(-^) 4 - (-^ f)dx) (69 

Vtot 2 Jl) Aoo 2 X A h iX 

where 

Aoo aa; = (Jfi-jV* = Aaa! 

1 2J* 2 
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and 

A h ix = (^!^)V3 = Dl 

1 2Vtot 1 

The value of A M 2 x has the sense of the spectrum width when the cross influence and the exhaustion of th 
heterogeneous centers of this very sort are neglected. The value of A h \X has the sense of the spectrum widt 
when the vapor consumption by the droplets is neglected. 

In addition one can easily prove that 

\ N 2 {i) - N 2 {j) |>0 



dx 



and 

' N 2{i) -N 2{j) \<0 



di]tot 1 

for i, j > 2. Hence, it is easy to show that 



N 2 (2 ) - N 2 
No 



< 0.015 



by the calculation of N 2 ( 2 ) (00) and N 2 (3) (00) at r\ to t 1 = 

The analytical approximation valid for the clear interpretation can be obtained if we notice that 



6 2 (2) (oo) = exp[-/. 2 A oc 2^(|(1 + (f^) 4 )" 1/4 + f (1 + (^) 3 )- 1/3 )] (70 

Vtot 2 2 I\ h ix 2 l\ h ix 

B = I exp(-x 3 )dy 
Jo 



with the relative error less than 0.02. 

The spectrum of sizes of the droplets formed on the centers of the second sort is the following one 



IW* 2 Z 4 T 2 TJtotl 3n , , "oo f Z , , X 4 X 3 

■ — r) ex P(~* z )exp(-/* 2 / exp(-(- ) - (- r) 



5.3.2 Situation hi > 1, ft 2 < 1 



The description of the process of formation of the droplets on the centers of the first sort can't be simplified. I 
has been already fulfilled in the previous sections. But the process of formation of the droplets on the center 
of the second sort is rather simple to describe. The supersaturation is absolutely determined by the vapo 
consumption by the droplets formed on the centers of the first sort. Hence, we have the following expressions 



6 2 = exp(-/* 2 — I cxp(-r 2 |^)dx) (72 

Vtot 2 Jo 



92 = U 2 J (* - exp(-T 2 ^-)9 2 dx (73 

The value of g 2 during the period of nucleation on the centers of the second sort can be estimated in th 
following way 

92 < =- 
1 2 

which is presented on the base of 

5 2 x <C A 2 x 

It is necessary to calculate only 9 2 . To calculate 9 2 one can get into account that the value of gi grows s 
rapidly that for the value J* exp(— T 2 ^-)dx one can show the following approximation 

[ exp(-r 2 fL )(& » 3 e(i - + ^« p( -r 2 |L Ml e ( ^ - 1) 
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j" exp (-r 2 £)dx « z6(i - ^) + 1) 

where z& is extracted by the condition 

t 2 

and is the Heavisaid's function. The last approximation solves the problem of the analytical calculation c 

5.3.3 Situation hi > 1, h 2 > 1 

Actually, the situation hi > 1, /12 > 1 has been already described in the previous subsection. Due to Aix <C 
one can't assume here that the inequality hi <C 1 ensures the pure exhaustion of the heterogeneous center 
without any vapor exhaustion and the consideration made in the previous section can't be simplified. 

5.3.4 Situation h x < 1, h 2 < 1 

From the first point of view it seems that the situation hi -C 1, /12 •C 1 has been already described in the sectio: 
5.2.4. We have to stress that the equation h <g 1 doesn't allow to state that the process of condensation i 
going at the constant value of the supersaturation. For the process of condensation on the centers of the firs 



sort we have previous expressions (62), (64). Analogous expressions (63), (65) for condensation on the center 
of the second sort can be violated. So, the process of condensation on the centers of the second sort can not b 
described on the base of the unexhausted value of the vapor supersaturation. 

The calculation of gi isn't necessary and only the calculation of #2 is essential. We have 

O2 (2) = exp(-/„ 2 — f exp((-(-^) 3 )dx) (74 
r\tot 2 Jo A h ix 

and the final value for 62 can be given by 

6 2 (2) (00) = exp[-/„ 2 A,, x x — —B] (75 

Vtot 2 

The monodisperse approximation is based on the evident chain of the inequalities 

Axi ~ Six <c Aix < Ax2 

or the monodisperse term isn't essential. 

We have investigated all possible situations.^] 

6 General approximate solution 

The reason why the general iteration procedure fails lies in the fact that we don't know the true expression fo 
gi(z). In the case of condensation on the centers of the separate sort we don't know this expression only whe: 
the spectrum is cut off by the exhaustion of the heterogeneous centers. So we don't know it in the situation the: 
the converging force of the operator in the expression for 6 is very strong. In the case of the several component 
the situation is different. We don't know every term in y\ 9j{ z )- Hence, we may come to the situation whe: 
according to the first iteration the spectrum is cut off by the exhaustion of the supersaturation initiated by th 
other component but in reality the heterogeneous centers of the other component are exhausted and and th 
droplets formed on these centers consume the vapor much more weaker (simply due to their quantity). So, it i 
necessary to invent the new more precise expression for gi which allows to calculate the next iteration in 9. 

6 It doesn't follow that all of them can be really reproduced in the nature. 
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6.1 Monodisperse approximation 

As it is stated in section 5.1 the length corresponding to the cut-off of the supersaturation which leads to th 
cut-off of the spectrum by the exhaustion of the substance is practically one and the same for all sorts of th 
droplets (all T t have approximately one and the same order). Let us see the droplets of what sizes play th 
leading role in this cut-off. Analyzing the subintergal expression in the equation for gi we realize that th 
subintegral expression connected with the variation of the supersaturation is a very sharp function of x. It i 
less than the function 

Sbei = Q(z-x){z-xf (76 

and greater than the function 

s ab = Q(z - x)(z - xf exp(- r ' ^ f * j ^) (77 

Let us introduce the approximation for this function. We must extract the region of the sizes of the droplet 
which are essential in the vapor consumption. This consumption is essential when 

x w Ax (78 

where Ax is the length of the cut-off by the supersaturation. Certainly, the region of the sizes of the droplet 
which are essential in the vapor consumption must have the sizes rather small in comparison with A^ becaus 
the successful iteration procedure in the homogeneous decay is based on the fact that the droplets formed at th 
almost ideal supersaturation determine the process of formation of the spectrum. For the differential halfwidt 
8i/2 we have the following expression 

S1/2X - (1 - (79 
The integral halfwidth A 1 / 2 x can be obtained from the corresponding equation 

x 4 

N ess x 3 = /* i— noo (80 
where N ess is the characteristic number of the droplets obtained as N ess = /* iAi^xn^ which gives 

A 1/2 x = jx (81 

and it practically coincides with Si/ 2 x. 

Figure 1 illustrates the behavior of the values s a b and Sbei as the functions of x. It can be seen that the 
practically coincide. This lies in the base of the applicability of the first iteration as a good approximation fo 
the precise solution. 

So, the subintegral function s is now decomposed into the essential part where 

A t x 

and the tail where 

AiX 

x > 

~ 4 

We shall neglect the tail and due to rather a small size of the essential region we shall use the monodispcrs 
approximation for the droplets formed in this region. As the result we obtain the approximation for g(x) 

g{2) = m^) z s (82 

where N(z/A) is the number of the droplets appeared from x = till x = z/4. 
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As far as the spectrum is cut off by the exhaustion of the supersaturation in a frontal (sharp) manner th 
value of gi is unessential before z = AiX as a small one and after the moment of the cut-off it is unessentu 
also as there is no formation of the droplets. 

So instead of the previous approximation we can use 

N t (A lX /A) 3 

gi{z) = z (oo 

rioo 

The process of the exhaustion of the heterogeneous centers makes the subintegral function more sharp an 
the monodisperse approximation becomes at AiX even better than in the pseudo homogeneous situation. Bu 
the exhaustion of the heterogeneous centers makes the coordinate of the cut-off of the supersaturation greate 
than AiX and the monodisperse approximation becomes even more better at the moment of the cut off c 
the imaginary supersaturation. Certainly, we must use N(Aix/A) calculated with account of the exhaustio 
of the heterogeneous centers (but at the coordinate, obtained without any account of the exhaustion of th 
heterogeneous centers). 

Figure 2 illustrates the form of the spectrum in the " monodisperse approximation" . The case of the pseud 
homogeneous condensation is considered. Here two curves are drawn: the spectrum in the monodispers 
approximation f apP r and the spectrum in the first iteration f\ which can be considered as a very precis 
approximation. The first iteration has the more sharp back side (the front side in z scale) than the " monodispers 
approximation. Nevertheless the deviation isn't so essential. It can be eliminated by the perturbation theory. 

The concluding remarks concern the fact that we can obtain N{A i x/ / i) by the solution of the equations fo 
the separate condensation process because we need the lowest length of the cut-off. This length is given withou 
any cross influence taking into account due to the frontal character of the back side of the spectrum. 



6.2 Final iterations 

The solution of the system of the condensation equations is given by the following procedure. 
At first we must solve the equations for the separate processes 

ft = /* i [ (z- xf expi-Tijg-Widx = Gi( gi ,0i) (84 
Jo <p * 

Oi = exp(-/* i—^- \ cxp(-r i -^-)dx) = Si(gi) (85 
Vtot i Jo ®* 

for every i. 

This solution can be obtained by the iteration procedure described in section 4. The remarkable fact is tha 
in the second approximation we can calculate the value of 9( 2 ){Aix/4) 

9 i{2) (A iX /4) = exp(-/* ( AtX,i cxp (-^^z 4 )dx) (86 

Vtot i Jo 4$, 

6 i{2) (A iX /A) = eM-f* > — (^r 1/4 C) (87 

Vtot i 4<P* 

where 

,1/4 

C= / cxp(-z 4 )dx w 0.25 (88 
Jo 

and iVj( 2 )(Aja;/4) has the following value 

N i{2) (A iX /A) = Vtot l (l - e i{2) {A lX /A)) (89 

So the approximations for gi are obtained now. We must fulfil these calculations for every sort of the heterc 
geneous centers. Considering these approximations as the initial ones it is necessary to do only one step of th 
iteration procedure to get the suitable results. They will be marked as the "final" ones. 
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Now we remove to the general iteration procedure. The only thing we have to do is to calculate 6i finalise 
and N fi na i j(oo). We have to obtain these values due to the iterations 



Vi ftnai{z) = exp[-/« 4 / cxp( J - z°)dx 

Vtot i Jo rioo'v* 



(90 



N fi na i i(oo) = ?7 tot i 1 - exp(-/* i (^ ) (91 

?7tot » noo$* 

All these equations are valid under the reasonable separation of the heterogeneous centers into the sort 
when the centers with one and the same height of the activation barrier are considered as one sort. Meanwhil 
it is obvious that if we split one sort into very many subsorts then we can formally attain a wrong result. 

Let us suppose that the number of the heterogeneous centers is so great that hi ^> 1, i.e. the fall of th 
supersaturation leads to the interruption of formation of the droplets. We shall split this sort into so man 
subsorts that for every subsort hij <C 1, i.e. we can see the exhaustion of the subsort in the separate process c 
condensation. Moreover, we can assume that this exhaustion finishes before z attains the quarter of the lengt 
of the spectrum initiated by the fall of the supersaturation. Really as far as AijX is proportional to Vi/tm 
(where r]ij tot is the total number of the centers of the given subsort) and Sijx doesn't depend on this quantit 
the required property is obvious. Then after the summation over all subsorts we come to the conclusion that th 
total number of the droplets formed up to Aix/4 coincides with the total number of the heterogeneous center; 
This conclusion isn't valid. 

The reason of the error is that the width of the spectrum is much more small than the width in the process c 
the separate formation. Evidently the monodisperse approximation doesn't work at such distances. Ordinary 
under the reasonable definition of the heterogeneous sorts all characteristic lengths are different. So, ordinary 
there is no such effect in this situation. 

Now we shall correct the theory. Note that the monodisperse approximation leads to some already define, 
functional form for g and for the supersaturation. We already know this form and now we need only to establis 
the parameters in these functional dependencies. 

We shall choose the unique length of the spectrum Ax. At this very moment this value is unknown, but A. 
satisfies the following inequality: 

Ax < Aix 

for every sort of the heterogeneous centers. Then for g we have 
which leads to 

JV i(2) (Ax/4) = ry, tot (l - 9 i{2) (Ax/4)) (93 

0i (2 ) (Ax/4) = exp[-/, ^ f exp(--^- V N, (Ax / 4))dx] (94 

Vtot i Jo 4 > *^oo 

We neglect here for simplicity the weak dependence 1^ = Y on the sort of the centers, while the strong dependenc 
/*j on the sort of centers is taken into consideration. After the substitution we get the system of the algcbrai 
equations for Ni(Ax/4) 

Ni{Ax/A) = rji tot (l - exp(-/* f**'* cxp(--^— V N 3 (Ax/A))dx)) (95 

Vtot i Jo n oo ^ 

Let us simplify the last system. To calculate the integral one can note that 



f x 

/ exp(— x 3 )dx w x x < 1/4 
Jo 



(96 
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So the trivial dependence of the r.h.s. on Nj disappears. We come to 



71 /\ 7* 

N^Ax/A) = Vl tot (l - exp(-/* i-^—f)) (97 

mot i 4 

On the other hand we can utilize the sense of Ax as the halfwidth of the spectrum due to the fall of th 
supersaturation: 

^5^(^/4)11 = 1 (98 

J 

After the substitution we can get the equation for Ax 

Ax 3 s-^ rioo Ax T 

— J>«(l-« p(-A J __))_ = l (99 

Equation ( p7[ ) expresses A r ^(Ax/4) through Arc and solves the problem. 

Equation (|9^) gives the expression for g and, thus, for the supersaturation as the function of time. Th 
number of the heterogeneous centers can be found by (|94|). The final values can be found by the correspondin 
previous formulas. 

We can estimate the relative error by the pseudo homogeneous case. Then we shall obtain that the relativ 
error can be roughly estimated^ as 

JV,(oo) - N fmal ((oo) | < \A-B\ ^ Q Q2 (10Q 



Ni(oo) ~ A 

The investigation of the process of the essential formation of the droplets is completed. 



7 Concluding remarks 

As the result of the previous sections we know the behavior of the supersaturation and the behavior of th 
number of the free heterogeneous centers during the period of the intensive formation of the droplets. It allow 
to get the main characteristic of the process of condensation - the total number of the droplets formed on th 
heterogeneous centers of the different nature. The final formulas allow the interesting quantitative physics 
analysis which is missed only due to the lack of the volume. The evolution of the already formed spectrum i 
much more simple. The high accuracy in the determination of the supersaturation isn't necessary now. Due t 
the proportionality of the velocity of growth to the supersaturation only the high relative accuracy is necessary 
The description of the further evolution is analogous to the investigation of the homogeneous decay made i: 

The description of the further periods can be given with the help of the direct application of the monodispers 
approximation as in the case of the homogeneous decay of the metastable phase. We have the simple differentic 
equation for some isolated hydrodynamic element 

-^2 = C = ^*-2^N fina n(oo) — (101 

a at * — ' TLao 
i 

which can be easily integrated as far as the r.h.s. doesn't depend on time. All terms in this equation are know 
from the previous analysis. So there are no problems in the description of the process until the coalescence 



7 Precise calculations for condensation of the separate sort give the relative error of the monedisperse approximation as 0.022 
in the total number of the droplets. The summation over all sorts can not increase the relative error 
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Figure 1 



Behavior of functions st, e i and s a b- 
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IFigure 2 



Form of the spectrum in the "monodisperse approximation". 
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